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We present a detailed analysis of glass transitions induced by pinning particles at random from 
an equilibrium configuration. We first develop a mean-field analysis based on the study of p-spin 
spherical disordered models and then obtain the three dimensional critical behavior by the Migdal- 
Kadanoff real space renormalization group method. We unveil the important physical differences 
with the case in which particles are pinned from a random (or very high temperature) configuration. 
We contrast the pinning particles approach to the ones based on biasing dynamical trajectories with 
respect to their activity and on coupling to equilibrium configurations. Finally, we discuss numerical 
and experimental tests. 

PACS numbers: 



I. INTRODUCTION 

A universal feature of super-cooled liquids is the huge growth of relaxation time when temperature is decreased 
below the melting point. In a restrict temperature window the relaxation time-scale increases by many orders of 
magnitude. Starting from the microscopic dynamic time-scale of simple liquids (picoseconds), the relaxation time 
increases so steeply that below a certain temperature, Tg, it exceeds the observation time scale (hours). Eventually, 
for T < Tg, it is not possible to equilibrate the system anymore: the sample is frozen in an amorphous solid called 
glassii^. 

Despite long lasting theoretical and experimental efforts, a definitive explanation for the origin of the dynamical 
slowing down in supercooled liquids is still lacking. The main theoretical question is what kind of critical behavior — 
if any — is responsible for the slow dynamics. Is it of thermodynamic nature, e.g. due to a Random First Order 
Transition (RFOT) at a temperature Tk < or to an avoided defect-mediated transition^, or is it instead 

purely dynamical^"—? In the case of usual second-order phase transitions, critical properties emerge and can be 
experimentally resolved within the critical region only, i. e. when time and length scales become much larger than the 
their microscopic counterparts. Unfortunately, in the case of the glass transition, all evidences point toward a modest 
increase of the length-scales associated to collective behavior, even close to T^. This is likely due to an exponential 
relation between time and length-scales different from the usual power law one found in standard critical phenomena^. 
Indeed, approaching Tg the relaxation time increases by several orders of magnitude, but typical length-scales increase 
by no more than a few inter-particle distances (ten at most). The consequences of this state of affairs is that: (1) it 
is not possible to observe genuine critical behavior, (2) testing and comparing different dynamic and static theories is 
very difHcult, (3) it is not known whether a true phase transition takes place at a finite temperature Tk below Tg or, 
instead, the transition is avoided or present at T = only. 

We recently proposed a way to partially solve these problems and to find and characterize the ideal glass transition at 
Tk , if present^. We proposed that ideal glass transitions can be induced in supercooled liquids at temperature higher 
than Tg by a suitable random perturbation: The idea is to block a fraction c of particles in the positions they have in 
a chosen (but arbitrary) equilibrium configuration at temperature T and study the thermodynamics of the remaining 
free particles. Using the RFOT approach we found that in a finite range of temperatures, T e [TK,Th], a true 
ideal glass transition, dubbed Random Pinning Glass Transition (RPGT), is expected to occur at finite concentration 
ck{T) of pinned particles. This transition for c t ck{T) has the same features as the ideal glass transition for T ], Tk- 
However, it can be studied much more thoroughly. The major advantage is that, as we shall explain later, the ideal 
glass phase is known in this case and, hence, it is possible to equilibrate the system, not only below ck (T) in the liquid 
phase, but also above ck{T) in the amorphous glass phase. For example, a standard local Monte Carlo or molecular 
dynamics can correctly probe the equilibrium measure on time scales not diverging with the system size in the glass 
phase. Therefore, for the first time, the ideal glass transition can be approached both from the liquid and from the 
glass phase. Showing the existence and the critical properties of the transition is not left to doubtful extrapolations 
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but can be ascertained by using the finite size scaling machinery developed for standard phase-transitions. Moreover, 
information on the presence, the position in temperature and the critical properties of the glass transition of the 
unperturbed system might also be numerically or even experimentally grasped by extrapolating at c = the analysis 
performed at finite c for T > Tg. 

The main aspects of the physical mechanism behind random pinning glass transitions and some results have been 
announced in Q. In this paper we present a comprehensive derivation and analysis of the theoretical predictions, 
including in particular the phase diagram in the c — T plane and the critical properties. The latter have been obtained 
through a non-pcrturbativc renormalization group method that we describe in detail. In this work we also compare 
our approach to alternative ones developed by Franz and ParisiiS and Chandler and Garrahan and co-workers2iiiii^. 
Despite similarities, especially concerning the phase diagrams, these two approaches turn out to be quite different: in 
particular they predict a first-order phase transition line and not a hona-fide glass transition line, as in our case. The 
difference between pinning from an equilibrated and from a random configuration is also addressed in this work. We 
shall show that the latter procedure leads to quite different phase diagram and physical behaviors. 

The procedure of pinning particles from an equilibrium configuration plays a key role in our work. It appeared in 
the literature already a decade ago in and After the suggestion in [l^, it became an important tool for 

testing theoretical predictions related to the ideal glass transition and probe medium range amorphous order— lii. A 
number of different geometries for the set of pinned particles have been studiedi^"— : the cavit y, t he two-wall and the 
one-wall geometry and pinning a fraction of particles chosen at random. Although already in islfl^ it was observed 
that this kind of procedures induces a significant slowing down of the dynamics of the remaining free particles, only 
recently it was claimed by numerical simulations that the impressive increase in relaxation time may be related to a 
true divergence^i^. In our previous work Q we show that this divergence is theoretically expected within RFOT 
theory and it is associated with an ideal glass transition, analogous to the one at Tk, occurring even at rather high 
temperatures, in some cases even close to the onset of glassy dynamics. The predictions of approaches other than 
RFOT theory concerning the effect of pinning particles appear to be quite different ranging from a gradual cross-over 
(instead of a singularity), as it is the case for some dynamic facilitation models^i, to no characteristic effect at all. 
For this and other reasons, the pinning particles procedure is a natural framework where different theories of the glass 
transition can be fruitfully challenged. 

The structure of the manuscript is the following: in Sec. 2 we review the phenomenological approach proposed in 
Ref . to explain random pinning glass transitions. The bulk of the paper that contains the mean-field analysis, the 
RG study and a toy model corresponds to Sees. 3, 4 and 5 respectively. Sec. 6 is devoted to a comparison between 
the random pinning approach and Franz-Parisi and Chandler-Garrahan ones. Finally, in Sec. 7 we present a general 
discussion, highlight the main advantages of the pinning particles procedure to study the glass transition and we 
suggest numerical and experimental tests. 



II. THE RANDOM PINNING GLASS TRANSITION: A PHENOMENOLOGICAL APPROACH 



In this section, we recall the phenomenological arguments we presented in [9| to explain random pinning glass 
transitions. Since our basic framework is RFOT theory we have first to present its main physical ingredient, which 
is the competition between two thermodynamic quantities: the configurational entropy Sc{T), which is the entropy 
per unit volume associated to the multiplicity of the amorphous phases in which a super-cooled liquid can freeze, 
and the interface free-energy cost AF of surfaces separating regions arranged in different amorphous phases. On 
length-scale i these two quantities respectively scale as Sc{T)£'^ and T£^ , where 9 is equal or less than d — 1 {d is 
the spatial dimension), T is the so-called amorphous surface tension. By lowering the temperature Sc{T) decreases; 
it is assumed to vanish with non-zero first derivative at a finite temperature Tk>^. This is at least what one finds 
in mean-field models^, approximate computations^^ and in extrapolations from numerical dat a^^'^'' . By comparing 
the scaling with £, one finds that on small length-scales the interface free-energy cost prevents the system from freely 
rearranging in different amorphous phases; the entropic gain due to the rearrangement instead becomes dominant on 
large length-scales. By balancing the two terms one finds a characteristic length £ps = (T/Tsc)^/*- called point- 
to-set. For £ < £ps, it is preferable for the system to be localized in a single given amorphous phase, characterized 
by a specific hidden amorphous order. For £ > £ps, the amorphous order is disrupted and the system is organized in 
a mosaic state, a kind of micro-phase separated state in which the number of competing phases is actually huge. The 
relaxation time is governed, in this picture, by cooperative activated rearrangements on the length-scale £ps- These 
involve barriers that scale with £ps and, hence, lead to a relaxation time r ~ TQexp{A£pg/T). Thanks to the power 
law divergence of £ps at Tk, the RFOT theory obtains a super- Arrhcnius law for the relaxation time-scale of the 
system and, for certain combination of the exponents 9 and ip, the Adam-Gibbs law. Note that there is not yet a firm 
microscopic derivation of the values of the exponents 9, ^. Several recent works have been devoted to this issue^^"— 
and heuristic arguments put forward in 0] suggest 9 = 4> = |- See the recent book[3l| for more details on RFOT. 
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Let's see how does the picture change when a fraction of particles is blocked from an equilibrium configuration. The 
exact procedure is the following: take an equilibrium configuration, pick at random a fraction c of particles, pin them, 
i.e. do not allow them to move anymore and then study the (static and dynamic) physical behavior of the remaining 
unconstrained particles. Of course both their configurational entropy density, Sc{c,T), and their amorphous surface 
tension, T{c,Y), are affected by pinning. However, we expect that the dominant change, at least close to Tk, is in 
Sc{c,T) since the configurational entropy is small close to Tk and the extra constraints due to the blocked particles, 
diminish it even further. Indeed, any blocked particle affects the possible configurations of its neighbors leading to a 
microscopic decrease of the effective configurational entropy density. For small c let us approximately write the new 
configurational entropy density as Sc{c,T) ~ Sc{T) — cY{T). Increasing c the configurational entropy monotonically 
decreases. For temperatures close to T^, Sc{c,T) is expected to vanish linearly as Sc{c,t) ~ Y{T){ck{T) — c) at the 
critical fraction ck{T) ~ Sc{T)/Y{T). 

By repeating the usual RFOT theory arguments but now with a configurational entropy that vanishes increasing c 
(instead of decreasing T), one obtains the physical picture of an entropy vanishing transition at ck(T), analogous to 
the one at Tk. In the pinned system, the effective competing terms are T{c,T)i^ and Sc{c,T)£'^ with T(c, T) ~ T(r) 
and Sc{c,T) < Sc(T). As a consequence, the liquid phase in presence of pinned particles is characterized by much 
larger cooperative rearranging regions of size £ps{c,T) = (T(c, r)/Tsc(c, T))^/'^'^"^' ^ £ps{T). The mosaic length- 
scale of the randomly pinned liquid grows when c increases and eventually diverges at ck (T) . Because of the relation 
discussed before between £ps and the relaxation time, one expects approaching to the entropy vanishing transition 
a generalized Vogel-Fulcher divergence r ~ Toexp{£psic,T)'^' /T) = tq exp(A(T)/(cif(T) - c)'''/^''-^)). At higher 
temperature this picture breaks down for the effective amorphous surface tension term T(c, T) changes. The reason 
is that fluctuations are so strong at high temperature that metastable states are not defined anymore, i.e. the 
interface- free energy cost T(c, T) is expected to vanish in the high temperature regime. Hence, random pinning glass 
transitions should only exist in a finite temperature range Tk < T < T^ where T^ is defined as the temperature such 
that lim^^y- T{ck{T),T) = 0. 

One of the main conclusion one can draw from the previous arguments is that by approaching random pinning glass 
transitions by increasing c toward ck{T) one should find the same physical behavior observed (or expected) for super- 
cooled liquids increasing T towards Tk- Instead, beyond the transition, there are important differences. In particular, 
whereas the ideal glass phase below Tk is completely out of reach and can be only found for small systems by solving 
a daunting optimization problem the one beyond ck is easy to obtain; it is given by the initial configuration used 
to pin particles (plus vibrations around it). The physical picture is that below ck the ideal glass is just one out of 
the very many states sampled by the liquid; at ck when the configurational entropy vanishes it is the only one to 
remain available. All other amorphous states have a higher free energy for c > ck{T). (This picture is quite similar 
to the one worked out for " quiet planting" in random constraint satisfaction optimization problcmsSS,) . A particularly 
instructive exam ple where this mechanism can be found at work is provided by the (pinned) Random Energy Model, 
see (and also j33l | ) . 

Note that the existence of a preferred state related to the initial configuration makes clear that pinning from an 
equilibrium configuration is very different from blocking the particle positions completely at random. In the latter 
case, there is no preferred configuration. In consequence, pinning acts like an external uncorrelated quenched disorder 
and finding the ideal glass state above ck is expected to be as difficult as finding it below Tk- We shall discuss in 
detail this point and the comparison between the two procedures later. 

III. MEAN FIELD THEORY OF RANDOM PINNING GLASS TRANSITIONS 

By now it is well established that the study of the disordered spherical p-spin completely connected model with p > 3 
provides a sketchy mean- field theory for the physics of supercooled liquidsii^i^ir— and RFOT. The mean field replica 
solution of this particular spin model displays an entropy vanishing transition at Tk, where the system undergoes 
a one-step Replica Symmetry Breaking (1-RSB) transition^. In this section we use this model to obtain a mean- 
field theory of RPGTs. Instead of pinning particles we block a fraction c of spins from an equilibrium configuration 
at temperature T. For comparison, we also work out the same analysis when the spins are blocked from a high 
temperature (essentially random) configuration, which mimics the case of particles blocked in completely random 
positions. 

In the context of combinatorial optimization, a similar analysis has been performed in Reflssl. 
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A. The spherical p-spin completely connected model 

Let us first recall definitions and notations for the usual spherical p-spin model (i.e. without pinning). Its Hamil- 
tonian reads: 

H = - ^ Jij...kSiSj . . . Sk (1) 

where the sum is performed over all the possible groups oi p spins in a system of N spins, Jij,,,k are i.i.d. Gaussian 
random variables with zero mean and variance aj — pl/2NP~^ , and the spin are continuous variables bounded by 
the following spherical constraint: N = X^il^i ^1- The mean intensive free-energy, T{T), averaged over the possible 
choices of the disorder, is obtained introducing n replicas of the system. Performing standard replica manipulations^^ 
we find (in the following we rescale T in units of fc^ and use the notation ~~ for the average over the quenched 
disorder) : 

F{T) = -^WZ = lim = - lim lim log / dQa,beM~NS{Qab;T)] (2) 

where Qab is a n x n symmetric matrix and 

1 " 1 
S{Qab) ^ ~^^Y.Qlb- l^JogdctiQab) . 

The saddle point solution for Qab in ([2]) is a matrix with n/m non-zero diagonal blocks of m x (to — 1) elements equal 
to qi, and m diagonal unitary elements. The saddle point equations determine the equilibrium value of qi and to. At 
a temperature Tk the free-energy has a singularity below which the equilibrium value of qi jumps from to a non-zero 
value gf"^ and to starts to decrease from 1, this is the signature of the IRSB transition. The other characteristic 
temperature in the model is the dynamic transition temperature T^. It marks the presence of a dynamical transition 
akin to the one found in Mode Coupling Theory; Td can be obtained as the highest temperature where one finds a 
stationary point of the action characterized by a non-zero value of qi for m = l^i^S. 

In the next two subsections we extend the usual replica computations in the case where a finite fraction c of spins are 
blocked according to the two different protocols cited above. 



B. Random pinning from equilibrium configurations 

The free-energy of a p-spin model with a finite fraction of spins blocked from an equilibrium configuration can be 
obtained through a double average process. For each choice of the quenched disorder, we have to average over the 
possible different equilibrium configurations used to pin the spins; then we have to average over the possible different 
choices of the couplings: 

W,c) = -^^P^(C/,T)logZ^^(c,T) , (3) 

Cf 

where P'' {Cf,T) = cxp{— 13 H'^ {Cf))/Z'^ {T) is the equilibrium probability of the configuration Cf in the free system at 
temperature T (= 1//3); Z^^{c,T) is the partition function of the system where cN spins at random are constrained 

to be in the same state than in the reference equilibrium configuration Cf, and Z"^ {T) = exp(— /3i7''(C)) is the 
partition function for the free system. Note that since the model is fully connected, there is no notion of space and 
the spin index can always be redefined in such a way that the the blocked spins corresponds the first cN ones. In 
the following we shall always follow this procedure. The average of the free-energy over the configuration Cf can be 
simplified introducing to replicas of the system, constrained to have cN identical spins, and sending to to 1 at the 
end: 

j:,^eM-PHJ{Cf)) ™™iTO-l°Sl^ Zc^cM-m^Cf)) )■ ^' 

The numerator in the RHS of (|1]) can be rewritten as a partition function, Z^^ , of to replica constrained to have the 
first cN spins identical: 

Z'^^^J2'^xp{-l3H'{Cf)){Zij-^-'= j2 oxp(-f2l3H-'{cA , (5) 

Cf Ci,...,C„ \ i=l I 
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where f2c and the prime above the sum stand for the constraint of having the first cN spins identical for the m rephcas. 
The average over the quenched disorder is performed introducing, as usual, n replicas and taking the n — > limit; 

- ^ hm - [log - log (Z^] . (6) 

Collecting the pieces together, we obtain: 



log (Z™ )" - log {z^Y 
Te(T,c) = — T lim lim lim . (7) 

Standard replica manipulations give the following result: 

{Z^^ J dQa^^b,eM-NS{Qa^b,;T,c)] (8) 



with 



/ n m 
\ a aaba 



S{Qa^b,;c,T) ^ ^{c+ {I - c)Q,^tJ'' + ^ QLb, ) -|(l-c)logdet(Q„„b, 



where the index Oq (and bp) runs over m different values for each of the n possible choices of a. Replica configurations 
identified by the same index a have the first cN spins identical, whereas replica corresponding to two different values 
of a have the first cN spins in uncorrelated configurations for T > Tk- The order parameter Qa^b^ is the overlap 
between the {l — c)N free spins of replica Cq, and 6/3 . The replica structure we focus on, and its physical interpretation, 
are identical to the one discussed by Monasson^: because of the pinning field induced by the pinned particles the 
m replicas inside each one of the n groups either fall into the same state, and are characterized by a high value of 
the intra-group overlap or fall into different states and are characterized by a small value of the overlap. Technically, 
this means that we consider a replica symmetric ansatz for the matrix Qa^bf, ■ Wc verified that taking a more general 
ansatz, e.g. IRSB, is not necessary as long as one focuses on temperatures higher than the glass transition temperature 
of the unpinned system. This is the regime we study in the following. Within the RS ansatz, Qa^bf, is a block matrix 
with n diagonal blocks having m diagonal unitary elements and m x (to — 1) off-diagonal qi elements. The off-diagonal 
blocks have only elements equal to qq. For this choice of the overlap matrix the action reads as follows: 

S{qi,qo,m;c,T) = -"x" [™("^ ^ (1 - c)(7i)p - m'^{n - l)q^] + 

-in(l - c) [(to - 1) log(l - qi) + (1 - i) log(l - qi + m{qi - qo))+ 

+ 7i log(l -qi+ m{qi - qo) + nmqo)] (9) 

and the saddle point equations over the action in the n — s- and to 1 limit give the following conditions on the 
parameter qi and qo: 

^(c+il-c)q,r-' = j^ (10) 

_P_„P-1 - 22 (-[-I] 

2T^% — (l-0n)2 y^^J 



solved by go = and a non-trivial qi. The vanishing of go is natural since it measures the overlap between two 
configurations corresponding to two systems with the same quenched disorder but different C/. In this case go is zero 
because two typical equilibrium configurations of the unconstrained system have zero overlap above Tk and, hence, 
the same is expected for the corresponding pinned configurations. 

At high temperature or low concentration c of pinned spins this set of equations is solved by a very small gi oc cP~^. 
However, for temperatures lower than Td{c), a solution characterized by a larger value of gi appears. This signals that 
below this temperature a metastable state lumping together configurations similar to the one used to pin particles 
starts to exist. Thus, the equilibrium dynamical evolution obtained starting from Cf is expected to show an ergodicity 
breaking. Since C/, as far as the unpinned spins are concerned, is statistically equivalent to any other equilibrium one, 
the dynamical behavior starting from any equilibrium configuration for the pinned system is expected to show just 
the same phenomenon. This allows us to conclude that at Td a dynamical Mode-Coupling-like transition takes place 
and that below Td the liquid state of the pinned system is broken in an exponential number of amorphous metastable 
states. Indeed, the value of Cd{T) obtained by solving the set of equations constituted by eq. (|10p and its derivative 



(12) 
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FIG. 1: Phase diagram for the spherical 3-spin model with spins pinned from an equilibrium configuration. The full 
line represents the thermodynamic transition ck{T) and the dotted lines correspond to dynamical transitions. The 

dynamic transition at higher c, Cd{T), is the usual Mode-Coupling transition line found in the study of the 
equilibrium dynamics. The other one, called Cdf(T), that we report for completeness can be only found through the 
study of non-equilibrium dynamics, as shown in^. Note that as discussed in the main text, and contrary to the case 
of completely random pinned spins, it is possible to continuously connect the liquid and the glassy phase through a 

closed path around the final critical point of the ck{T) line. 

coincides with the one found in the dynamical analysis presented in [4lj . 

The secondary minimum, which appears at Cd{T), becomes the global one when: 

-^(c-f(l-c)gif -(l-c)[log(l-gi) + gi]=0 . (13) 

This equation defines the line of thermodynamic transition ck{T). As noticed in the previous paragraph, since the 
metastable state associated to the initial configuration is statistically similar to all other ones, this means that at 
ck{T) the free energy of the liquid state becomes equal to the one of a given metastable state, i.e. the configurational 
entropy is zero. In consequence, the transition at ck is an entropy vanishing transition as the one reached for T — >■ Tk. 
For c > Ck the equilibrium phase of the pinned system is given by the state associated to the initial configuration. 
In this regime, there are still other states but they are characterized by a higher free energy. In App. |^ we present a 
more detailed general argument to explain why the equal value of the action at ck{T) implies Sc(c, T) = 0. In App. [C] 
we directly compute Sc{c,T) by generalizing the Franz and Parisi e-coupling approach to the case of a system with 
pinned particles and confirm that Sc(c, T) | for c t ck{T). 

The complete phase diagram for the p — 3-spin model is shown in Fig. [T] As shown, the RPGT line stops at a 
critical point at which the dynamical transitions merge with the static one. At this point the glass transition becomes 
continuous, i.e. the difference between the high and low overlap values vanishes. The two lines Cd{T) and Ck{T) have 
been obtained by a numerical solution of the corresponding two sets of equations. They had been already shown in 
the phase diagram of and [4l| . 

C. Random pinning from completely random configurations 

As already stressed before, it is natural to wonder how much the results presented in the previous section depend on 
the pinning protocol, in particular what is the difference between pinning from a random versus from an equilibrium 
configuration. Our expectation, based on heuristic arguments, was that the difference is substantial^. In this section 
we confirm this by repeating the analysis performed above when cN spins are blocked in a completely random 
configuration. 

In this case the average over the reference configuration C / is performed over a flat distribution. The computation in 
this case requires the simple introduction of n replicas all constrained to have the same (random) spin values on cN 
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FIG. 2: Phase diagram for a spherical p-spin {p = 3) model with spins pinned from a random configuration. The full 
lines represent the thermodynamic transitions and the dotted line marks the dynamic transition. There is a full low 
temperature region in the phase diagram characterized by a IRSB glass phase. This phase is separated from the 
liquid phase by a IRSB discontinuous-transition line ranging from (c = 0,Tr-) to a critical point (co,To) and by a 
IRSB continuous-transition line ranging from this point to zero temperature and c = 0.5. 



sites randomly chosen (we used the same notation than in the previous section): 



Tr{T, c) = -r lim lim °^ , (14) 



where 



and 



dQa^b, exp[-iV5(Q,„b, ; T, c)] (15) 



^«(Qa„^; c, r) = ^ 5] (c + (1 - c)Q,^^;f -2^^-^) logdet(g„„,,) . 

a/3 a^bjj 

This case of completely random pinned spins is equivalent to a p-spin model in which on each site there is a random 
field that is very strong and equally probable in sign with probability c and that is zero with probability 1 — c. As 
for the p-spin model in a uniform field, it is necessary to consider a IRSB ansat^, where Qa^b^ presents a block 
structure with n/m diagonal blocks where the non- unitary values are the off-diagonal qis and the elements in the 
off-diagonal blocks are all equal to go (we have verified that the IRSB ansatz is stable and there is no need for further 
RSB). By plugging this specific form of Qa^bf, into the action we find; 

5(gi, go, m- c, T) = -^n [(m - l)(c + (1 - c)q^Y + ^ - l) (c + (1 - c)q^Y] + 

^\n{l - c) [(1 - i^) log(l - q,) + - i) log(l - 91 + m(<zi - q„)) + 

+ j-^log{l- qi+m{qi- qo)+nqo)] . (16) 
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The saddle point equations for the parameters qi, qo, and m of the equilibrium matrix are the following: 



^{c+{l-c)q,r-'^- 
Zl m 



1 



1 



mqo 



1 - 91 1 - 91 + m{qi - qo) [I - qi + m{qi - qo)Y _ 

—2 c + 1 - c)qo)P = ■ 

2T^ (I - qi + m{qi - qo)y 



^ [(c + (1 - c)q^r - (c + (1 - c)qor] + 



1 log ( 

m VI - gi + m(gi - (Jo) 



1 - gi + m(gi - go) (1 - 91 + "lift - 9o))' 



(17) 
(18) 

(19) 



Since this problem shares many similarities with the problem of a p-spin model in a field, it is natural to find, as we 
do, a non-zero solution for qo and the presence in the phase diagram. Fig. [5J of a continuous transition line between 
the paramagnetic and the IRSB phased. This line corresponds to the points in the c — T plane where the previous 
set of equations admit a solution with qi — qo ^ 0. To obtain the continuous transition line, following the procedure 
of (36j . we developed the full set of equations for small qi — qo and we obtained the following conditions on qo, m and 
c for any T: 



2T' 



P{c+{l-c)qoY-^ 



go 



2T2 



llil-c){c+{l-c)qor-'^j^ 



(20) 
(21) 
(22) 



By solving these equations in the range where < m < 1 (corresponding to {p — 2)/(3p— 4:) < qo < 1) one finds the 
continuous transition line cj!(r). Out of this range, the transition becomes discontinuous in qi — qo and it is always 
characterized by the value m — 1. The equations p7|) . (fT8| . and ([T9)) . in the m — )• 1 limit, define the discontinuous 
transition line c^j^{T). In the same limit, also the dynamic transition line c5(T) can be obtained requiring that the 
derivative of (|17p with respect to qi is zero instead of satisying equation (jl9|) . All results are summarized in the phase 
diagram of Fig. [2] 



D. Discussion on the phase diagrams and comparison between the two pinning procedures 

We now present a contrastive description of the phase diagrams (see Figs. [1] and [2]) obtained by the two pinning 
procedures studied above. Let us first focus on the thermodynamic transition lines. In the case of spins pinned from 
an equilibrium configuration, we obtained a line of thermodynamic glass transitions which separates a liquid phase 
from an ideal glass phase. These transitions are associated to the vanishing of the configurational entropy in the 
pinned system. Two important and related features of these transitions are: (1) the transition line stops in a critical 
point, (co,ro), where it becomes a continuous transition (2) the glass phase can be obtained from the liquid phase 
without crossing any thermodynamic singularity. This is due to the particular kind of low temperature phase of the 
system, which corresponds to configurations close to the reference amorphous state and not to the usual non-trivial 
1-RSB ideal glass phase. Indeed, the overlap distribution of equilibrium configurations is trivial and characterized by 
a single peak in the glass phase in contrast to what happens in 1-RSB phases where it is bimodal, as we shall discuss in 
more detail later. Even if there is no periodic order in this phase, a non-evident amorphous order is directly imposed 
by the pinned spins (or particles). As a consequence, to make the difference with usual 1-RSB glasses, we refer to 
it using the evocative name "amorphous crystal" [(AC) in the notation of Fig. In the case of spins pinned from 
a random configuration, the glass phase is instead 1-RSB and it cannot be reached from the paramagnetic (liquid) 
phase without crossing, either a discontinuous 1-RSB transition phase, or a continuous one. This can be understood 
easily since the completely random pinned spins act as an uncorrelated quenched disorder (actually a random field) 
and, hence, the usual RSB glass phases [(G) in the notation of FigIS], naturally emerges at low temperature^. As 
shown in Fig. [3] where we compare the two phase diagrams, the discontinuous glass transition line emanating from 
c = 0, T = Tk is quantitatively similar for the two procedures for small values of c. Instead, after the critical point 
Co, To delimiting the final point of the discontinuous transition lines, the behavior is quite different since in the former 
case the line stops whereas in the latter the line continues until reaching zero temperature for much larger values of c. 
We now focus on the dynamical transitions. In the case of the first pinning procedure one finds (see Fig. [T]) a Mode 
Coupling Transition line, Cd{T), and also a new dynamic transition line Cdf(T). The latter can only be found in 
the out of equilibrium dynamics: the relaxation time to equilibrium after quenches from high temperature diverges 
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FIG. 3: Zoom of the rectangle enclosed by the dashed line in Fig. [5] and comparison between the thermodynamic 
glass transitions obtained using the two different protocols to block spins. Top panel: the full line refers to the glass 
transition induced by pinning spins from an equilibrium configuration; it stops at the critical point {co,To) at high 
temperature. Bottom panel: the full line refers to the glass transition induced by pinning spins from a completely 

random configuration. It is a line of 1-RSB transitions which arc discontinuous on the left of the critical point 
(co, To) and continuous on the right of it. Below the full line the system is in a one-step replica symmetry broken 

phase. 

approaching the c^j. (T) line from the right even though the decorrelation time of equilibrium correlation functions 
stays finite, see |4l| . In conclusion, quenches from random configuration lead to aging in the region between the two 
dynamic transition lines. The situation in the case of spins pinned from a random configuration (Fig. [5] and bottom 
panel of Fig. [Sj is quite different. The Mode Coupling Transition line, Cd(T), is present also in this case. After the 
terminating point co,To, the line coincides with the static continuous one and corresponds to continuous equilibrium 
dynamical transitions. Note that the glassy region in which systems quenched from very high temperature remain 
forever out of equilibrium is much wider in this case compared to the previous one. 

IV. RENORMALIZATION GROUP ANALYSIS OF RANDOM PINNING GLASS TRANSITIONS 

In order to go beyond the mean field analysis and study the critical properties of the glass transitions on the ck (T) 
line, we follow the Renormalization Group (RG) procedure proposed in [42| and extend it to the case of systems with 
pinned particles. This analysis allows one to obtain the divergence of length and time scales that cannot be obtained 
within mean field theory. 

A. Field theory description of a replicated liquid 

The starting point of the RG study is a replica field theory based on a Ginzburg-Landau action, which combines 
the one found within mean-field theory with a square gradient term that disfavors spatial fiuctuations of the overlap. 
It can be obtained studying the thermodynamics of m — 1 copies (or replicas) of the same liquid system coupled with 
a small attractive interaction whose amplitude is set to zero after taking the thermodynamic limi t^^'^° . It gives a real 
space description of structural liquids in terms of the usual order parameter for the glass transition: the overlap ^ 

between an equilibrium configuration C and a reference configuration C, also at equilibrium. In order to construct 
a RG theory of structural liquids, one has to consider local fluctuations of the order parameter by introducing a 
local version of the overlap, g-(r) (obtained through a suitable coarse graining over length scales A larger than 
the inter-particle distance, ocr^, and smaller than the non-trivial cooperative length scales of the system). In the 
replica language this quantity is represented by the matrix of local overlaps qabir) between different replicas a and b 
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of the system. The simplest Ginzburg-Landau functional obtained retaining the leading terms of the overlap in the 
field theory and belonging to the universality class of systems showing Random First Order Transition at is the 
following^: 

^ T E {dqa,{x)f + V{c,{x))\ (23) 



where 



/in u + w ^ y 4 \ " 

^ ^ 2^ [ —Qab + -^Qab ) - ^ qabqbcQca ■ (24) 



a.b—1 a,fc,c— 1 

The temperature dependence in the action is only introduced in the parameter t = '^7^° , where Tq is a suitable 
scale of the temperature. The other parameters of the potential are positive and independent from the temperature. 
Compared to the action obtained within mean-field theory, a new space dependent term appears in the action to 
account for the cost of spatial fluctuations of the order parameter. This real replica description of supercooled liquids 
allows one to obtain the properties of the metastable states from the knowledge of the replica partition function, 

Z{m) = J YlVqabix)cM~S[q]) . (25) 

ab 

The mean free energy of a typical equilibrium state and the corresponding configurational entropy respectively read^ 
/3F = —9 log Z{m) I dm and Sc = — m^9(log Z{m) /m) / dm. The number m of replicas should be analytically continued 
to 1 in the equilibrium liquid phase and to a value less than one in the ideal glass phased. At the mean-field level, i.e. 
by looking for the uniform saddle-points of (j25p . one finds that the order parameter qab is zero above a temperature 
Td, such that id = -j^, and that below another uniform solution appears with a replica symmetric (RS) structure 
Qab = Qea > for a b. By explicitly using the RS structure of qab one finds that qea is the secondary local minimum 
of V{q) = and that the configurational entropy per unit volume is V{qEA)- Finally, at a temperature Tx 

m—l 

such that tfc = there is a random first-order transition with a coexistence between a zero-overlap phase and a 
high-overlap one, zero latent heat and vanishing configurational entropy density. Below Tr-, the system is in an ideal 
glass phase characterized by a nonzero overlap matrix and a value of m less than 1. We have chosen, without much 
loss of generality, a simpler form of V{q)^ which corresponds to impose qea = 1 at Tk. In this case one finds that 

where by construction and 7V|^^g respectively are the mean-field dynamical and Kauzmann transition tem- 

peratures. The sub-index c = is introduced for later convenience to recall that these are characteristic temperatures 
of the unconstrained super-cooled liquid. 

The effect of freezing a fraction c of particles can be schematically included in the field theory by forcing the overlap 
between replica to be equal to qsA at a set of points associated with frozen particles. In reality, the effect is more 
complicated than that but this is irrelevant as far as the large lengthscale properties are concerned. Thus, we impose 
to the measure in (j25p the constraints qabix) = qsA in a random set of Poisson-distributed points characterized by a 
density cp, where p is the particle density. 



B. Migdal-Kadanoff renormalization scheme and mapping to the Random Field Ising Model 

As discussed in [43 | in order to go beyond the mean-field analysis of RFOT one has to use a non-pcrturbative 
renormalization group approach. Following (4^ we first discretize the field theory and then apply a Migdal-Kadanoff 
renormalization scheme. This RG method becomes exact on hierarchical diamond-like lattices. Such lattices are 
built iteratively by replacing each bond between sites by a fixed number of new bonds which, to mimic Euclidean 
d-dimensional lattices, is taken equal to 2'^. By construction one can then integrate out iteratively degrees of freedom 
in the inner bonds and obtain a renormalized effective interaction, see Fig. |4l After n iterations, the volume of 
the system, which is proportional to the total number of original bonds, is renormalized by a factor 2"'' whereas 
the "distance" on the renormalized lattice corresponds to 2" original bonds: this naturally fixes the length-scale 
after n iterations as = 2". In the case without pinned particles at the n—th step of renormalization all links are 
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FIG. 4: Elementary step illustrating the RG procedure on a hierachical lattice. The notation is explained in the 

main text. 



characterized by the same renormalized effective pair interaction, W^{q^ , q^), which satisfies a closed equation written 
in terms of the pair interactions, W;Jij(q^,q^) and W^":i{q^ ,q^), of the links connecting those sites: 



(27) 



ji,,ij = l,2'i-i a,b 



where the labels 1 and 2 denote the value of two renormalized sites from which emanate 2"^'^"^) original bonds, and 
the labels L and R indicate the j th left and right link respectively, see Fig. |4l At the nth iteration, the original lattice 
is replaced by a renormalized one where the unit length is £n and the pair interaction between sites is W^(q^,q^). 
Freezing particles at random does not change the elementary RG equation (j27p but it requires that in the first 
iterative equations corresponding to n = 1, qab is fixed equal to qea with probability c for each given intermediate 
site^. As discussed in the following, the presence of these random constraints considerably complicates the analysis 
because the W„'s become random variables. In order to obtain a tractable problem we use the insight gained in 
the work [42|. Without frozen particles we have found that the nature of the RG flow is similar to that found for 
first-order discontinuity fixed points in which, generically, there are two essential couplings: the field h favoring one 
phase with respect to the other and the coupling J opposing spatial variations of the order parameter. In our context, 
h and J respectively correspond to the configurational entropy Sc favoring the zero overlap phase and the interface 
free-energy loss Y between a high-overlap and low-overlap phase. Guided by this result, we repeat in the following 
the analysis performed in [42j . "projecting" the functional MKRG equation on simpler ones, where the matrix qab 

is a two-state variable, qab = or qab = qsA for all pairs a,b. As in we focus on (91,52) = 



m— 1 

and V{q) = , and we approximate the full integral in (P7)) by a steepest-descent calculation (see [42| for a 

m— 1 

discussion of this approximation) . The iteration equation simplifies to 

W:,{quq2)^ J2 |w?f„i(qi,g) + W?,i!!i(<Z,g2) + ng)|, (28) 

, , q=0,qEA ( ) 
jt,n=1.2<'-i 

for qi and 92 that take the values and qEA only. The initial condition is W^Q^'"(qi, (72) = (91 — 92)^/2. We have 
checked that the results of [42j continue to hold within this approximation, as required. Because of the frozen particles, 
the above equation is modified at the first stage of iteration, n = 1 . For a fraction c of randomly chosen initial sites^, 
instead of minimizing over q, one just evaluates the expression within brackets in (|28p for q ~ qEA- In consequence, 
even though the RG equations do not contain quenched disorder, their initial conditions do. Thus, the W^{qi,q2)s 
become random variables and equation (|28p defines the flow equation for their probability distributions. Actually, it is 
useful to use the terminology of spin systems and interpret eq. (j28p as a minimization equation for a local contribution 
to the global Hamiltonian. Since a generic pair interaction W{qi, 52) can only take four different values depending on 
the values of 91,(72 one can rewrite it, without loss of generality, as 

2W{qi,q2) = - Ji2^1^2 - hiSi - h2S2 + C Si= ±1, 52 = ±1 

where up (down) spins correspond to high (low) local overlaps and C is a constant (the two in the RHS has been 
added for consistency with the notation used in ji^). Using this dictionary eq. (|28l) becomes identical to the zero- 
temperature MK-RG equations for the random-fleld Ising model (RFIM) in an external negative magnetic field equal 

to ~2V{qEA)- The initial condition for the MK-RG equation, VFq^"" (91, (72) = {qi — <Z2)^/2, leads to the initial 

2 

condition for the flow: — — 0,h2 — 0. At the first iteration, n = 1, for a fraction c of randomly chosen 

initial sites, instead of minimizing, one just evaluates the expression within brackets in (j28|) for S = I. This introduces 
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quenched randomness in the iterative RG equation: it is equivalent to apply a very large magnetic field, much larger 
than the ferromagnetic coupling on randomly chosen sites. In consequence, we find that the MK-RG equations for 
the replica field theory maps on the zero-temperature MK-RG equation for a RFIM with an unusual random-field, 
whose distribution is bimodal and asymmetric. It consists in a peak of weight 1 — c centered on —2V{qEA) and a 
peak of weight c centered at a positive value large enough to fix the spins in the up stated. In this scheme, the glass 
transition corresponds to the first-order phase transition between the negative (low overlap) and the positive (high 
overlap) magnetization state of the RFIM. The transition is controlled by the uniform negative field, which is equal 
to minus the configurational entropy, and the fraction c of frozen up spins. For each value of the negative magnetic 
field a transition is induced by changing the number of frozen up spins. Clearly, the larger is the absolute value of the 
uniform magnetic field, the larger is the critical value of c needed to counterbalance it. It is only when the fraction 
of frozen spins is zero, c = 0, that the transition takes place at zero magnetic field (i.e. —2V(qEA) = 0). 




(a) 1st iteration of RG flow (b) 3rd iteration of RG flow 




(c) 4th iteration of RG flow (d) 13th iteration of RG flow 



FIG. 5: Probability distribution of the coupling constants J and h, P{J, h^, hji)\f^^^j^^^f^, along the RG flow in the 
case Tdlj.^Q / TkI^^q ~ 1.4 at T = 1.35 7V|c=o ^^"-^ ~ 0.207435. The distribution, initially being a pair of delta 
functions one with weight 1 — c at small positive field and the other with weight c at a large negative field, is 
constituted by five delta functions at iteration one (see Panel [5a|) where only one delta function corresponds to small 
positive field (note that its weight is still the largest one) and by a smooth distribution formed by two almost equal 
parts corresponding to positive and negative fields for iterations larger than three (iteration three in Panel [5bl 
iteration four in Panel [5c|) . Eventually one of the two parts shrinks and disappears at £ = £ps = 2" (iteration 

thirteen in Panel I5d[) . 



C. Renormalization Group flow and phase diagram beyond the Mean Field approximation 

In the literature there are already several analyses of the MK-RG eqs. for the RFIM. However, the particular 
random field distribution that results from the mapping discussed in the previous section was not considered in any 
of those studies and, hence, required a new analysis that we present in the following. We have numerically studied 
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the flow of the joint probability laws of J, and hn by solving the iterative equations by population dynamics^. 
The initial conditions for the RG flow for a given link are J = and h]^ = hn = 0, as discussed in the previ- 

ous section. At the first step of the RG flow a fraction c of spins are frozen in the up state. This is equivalent to 
introducing random strong positive local fields that, if present, polarize the central spins in the up position. For any 
central site, the field is added with a Poisson probability characterized by a first moment equal to c. Starting from 
this random initial condition one then iterates the MK-RG eqs. We computed renormalized couplings and fields by 
using populations of 10^ variablesi^. 

Compared to the case studied in [i^] the analysis is complicated by the fact that one has to follow the flow of proba- 
bility distributions and not just of few coupling constants. Fortunately, even though initially these distributions are 
bimodal with one very sharp peak at the first steps of the RG flow, soon enough they become smooth and broad 
and they can be well characterized by their mean and variance, see Fig. [SJ The glass transition, like in (42j . is found 
studying when the high overlap-positive magnetization state becomes stable. We find that the values of the mean 
and variance of the couplings and the fields start growing along the renormalization group flow when the temperature 
is low enough. In the liquid phase there is always a characteristic scale above which replica between different sites 
decouple because the average coupling eventually becomes smaller than the average field, as shown in Fig. [7] (top 
panels) where we plot the ratio h/{h + J). The iteration n* at which this happens identifies the characteristic mosaic 
length-scale £ps — 2" , called point-to-set in the glass literature^. The longer replica remain coupled, the closer 
to the glass transition the system is and the larger £ps is. The line of random pinning glass transitions is obtained 
studying the critical concentration ck{T) at which £ps diverges for any given T. In Fig. [S] we plot the phase dia- 

gram obtained for y^ "" ~ 1.4. We recall that 7d|^^Q and Tk|^^q respectively are the mean-field dynamical and 

Kauzmann transition temperatures for the unconstrained system obtained from the analysis of V{q), i.e. in absence 
of fiuctuations (see eq. [26|) . Note that, as discussed in [i^l, Tr-|c=o '^'^^ renormalized by fluctuations because of the 
saddle-point approximation, whereas does not exist anymore beyond mean-field theory. By comparing Fig. [1] 

and Fig. [S] one finds that the mean-field phase diagram is drastically modified by taking into account fiuctuations. In 
particular, the endpoint of the RPGTs line, T^, turns out to be smaller than the mean- field mode-coupling transition 
temperature in absence of pinning, T^I^^q, contrary to what happens within mean-field theory, see Fig. [TJ Finally, 

the analysis of liquids with different fragility show that the ratio -pp-r — increases when „ decreases, i.e. for more 
fragile liquids. 

The continuous line in Fig. |6] is obtained using the relation ck(T) cx Sc(T), which was motivated by phenomcno- 
logical arguments at the beginning of this manuscript. It works surprisingly well. The RPGT line ends in a point 
corresponding to the concentration and the temperature c^, T^. Within the magnetic analogy, the physical reason for 
this phenomenon is that by increasing T the fraction of random spins frozen in the up position, which is needed to 
counterbalance the large negative uniform magnetic field (the configurational entropy), increases too. The net effect 
is an increase of the quenched disorder. As it is known, a strong enough disorder destabilizes long range order in 
the RFIM even at zero temperature. The endpoint of the line of RPGTs is precisely related to this effect: above 
Th there is no transition between high and low overlap phase because the quenched disorder introduced by random 
pinning is so strong that long-range order is wiped out. In two dimensions, as it is well known for the RFIM, even an 
infinitesimal disorder destroys long-range order. For this reasons no RPGTs are expected in dimensions lower than 
three. 

In the following we discuss thoroughly the critical properties of RPGTs in three dimensional systems. We separate 
two cases depending on the closeness to the endpoint Cfi,Ti^. Far from it the system behavior is controlled by the 
discontinuity fixed point as already found in [4^ without pinning, i.e. the critical properties at RPGTs are the same 
of usual glass transitions (without pinning). Instead, close to the endpoint the zero temperature RFIM fixed point 
intervenes and controls the critical behavior on time and length-scales that diverges approaching Ch^Th. 

1. Critical regime I: low pinning 

Let us focus first on temperatures slightly above Tk. In this case, by increasing c toward its critical values ck{T) 
we find results very similar to the ones obtained decreasing T towards Tk at c = 0. 

The average value of the coupling at the n-th iteration, J^, increases as y(c, T)£'^~-^, where Y{c, T) depend on T and c 
but remains of order one in the interesting ranges of these two parameters. On the other hand, the average value of the 
field increases as Sc{c,T)i''', where Sc{c,T) is by definition the configurational entropy (or the effective magnetic field 
in spins language) coarse-grained over lengths larger than the distance between pinned degrees of freedom. After some 
RG iterations the variances of the field and of the coupling become much smaller than the corresponding averages, see 
lower left panel of Fig. [71 This means that disorder is irrelevant on large length-scales and, hence, one finds back the 
same physical behavior discussed in [Z^l for unconstrained super-cooled liquids. In particular, the scalings Y{c, T)l'^~^ 
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FIG. 6: Phase diagram obtained by the renormalization group approach for 



1.4. The X-axis reports the 



amount of the concentration of pinned particles in the system. The squares are the results obtained by RG for the 
ideal glass transition line Ck{T), or Tk(c). The continuous line corresponds to the phenomenological result c/^ oc Sc. 
The endpoint of the glass transition line takes place at Ch,Th. Note that the temperature T/j is smaller than the 
mean- field Mode Coupling Transition temperature Tdj^^p for the unperturbed liquid (found using the Ginzburg 
Landau action with Td\^^Q / TkI^^q — 1.4) and, therefore, of the value T^^^ > T^I^^q found within mean-field 

theory. 

and Sc(c, r)£'* hold until one reaches values of £ such that the couplings become of the same order of the fields. Above 
this characteristic length-scale the average field starts to dominate and the average coupling decreases (it eventually 
vanishes for i — oo). Physically this means that the effective degrees of freedom of the renormalized system become 
non-interacting. This cross-over identifies the point-to-set length scale, ips, which physically corresponds to the linear 
size of regions able to decorrelate from their boundary^. Our numerical analysis shows that Sc (c, T) vanishes linearly 
by increasing c at fixed temperature, i.e. Sc(c, T) cx {ck{T) — c). In consequence the scaling of the point-to-set length 
is: 



£ps oc 



ck{T) 



As already discussed in [42|, on the length-scale £ps the renormalized system is like a liquid at its "onset temperature" 
where the (renormalized) point-to-set length is equal to one. The difference with a normal liquid are the values of 
the renormalized surface tension and configurational entropy, Y{c)i'''~^ and Sc^''- These are very large compared to 
T, i.e. the liquid is at very low T compared to the typical scale of the interaction. This suggests that the relaxation 
time could be obtained by assuming an Arrhenius law at the scale £ps and leads to a generalized Vogel-Fulcher law: 



log T oc 



Yi 



d-1 

PS 



1 



T 



icK{T)-c) 



d-l 



This critical behavior is valid for all temperatures above Tr- and below Th when approaching the RPGT by increasing 
c. However, its region of validity shrinks approaching because another fixed point starts to dominate the RG fiow. 
Discussing in detail the increase of time and length scale close to the endpoint is the aim of the following section. 



2. Critical regime II: strong pinning and the zero RFIM temperature fixed point 



As previously discussed and shown in the phase diagram, the fraction of pinned particles needed to induce a RPGT 
increases with temperature. This lead to an enhancement of the quenched disordered and, hence, of the variance of 
the random field in spins language. Previous works on the RFIM with Gaussian and symmetric bimodal disorder 
distributions have shown^=^ that, even at zero temperature, long-range order, and therefore the first-order transition 
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FIG. 7: RG flow of the coupling constants approaching the transition, i.e. for ck{T) — c — > 0^, projected on the 
planes a = and h = 0. The left panels correspond to a temperature very close to T/f lc=o' central ones to a 
temperature intermediate between Tk\^^q and T^, T/ Tk\^^q — 1.35. The right panels correspond to a temperature 
very close to T/j. The reported data, h, J, arc the average values of /i'^, J'^ and a is the variance of hl^. In the three 
cases, the projection on the ct = plane leads to the same flows obtained in^: J increases until the point to set 
length is reached. At that point h and J are of the same order. On larger length-scales J decreases rapidly and h 
continues to increase. In the left panels, the projection of the flow on the h = plane shows that the value of a 
remains small: in this regime the disorder does not change the sealing of J, h which are the ones typical of usual first 
order transition. For the central panels, a always decreases rapidly approaching a decreasing limiting curve for 
ck{T) — c — > 0^. For ^ :2> ^, when the renormalized a is small enough, the flow recovers the properties due to 
discontinuity fixed point, see left panels. For the right panels, the flow is dominated by the RFIM critical point. In 
this case a first increases as J because the flow is attracted by the zero temperature RFIM fixed point. The limiting 
curve, for ck{T) — c — > 0^ is a constant. This holds until the point to set length, £ps is reached. 



between negative and positive magnetization, are wiped out by a large enough disorder. This is indeed what we find 
(see the phase diagram): above a certain temperature, T/i, the disorder fluctuations induced by freezing particles are 
strong enough to make the glass transition disappears. The disappearance of long-range order due to the quenched 
randomness is controlled by a zero temperature flxed point different from the discontinuity one we studied in the 
previous section. Both are expected to play an important role close to Th, as we shall indeed show in the following. 
This connection with the RFIM could actually be anticipated on general grounds and has also been found in a 
complementary way based on a one-loop RG treatmentil. 

We now present our numerical analysis of the RG flow close to the endpoint of the glass transition line. Note that 
within the mapping to the zero temperature RFIM, the transition at Th is analogous to the one taking place at Ahc, 
the critical value of the variance of the disorder above which no long-range order is possible even at zero temperature. 
As a consequence, as it happens for the RFIM we expect the glass transition to become second order at T/j and 
CK{Th). Right at this point, the scaling properties are completely dominated by the zero temperature RFIM fixed 
point. The most important effect is that the average and the variance of both fields and couplings remain of the same 
order—, i.e. disorder remains relevant on all length-scales. We find that the average values and the variances of j', ft,' 
ah scale as with 9 = 1.491(4) for d = 3. This is the usual value found by MK-RG for the RFIM^, to be compared 
to the one obtained by simulations 6 = 1.49(3}^. When T is close, but smaller than T^, the critical properties of the 
glass transition become quite complicated since the RG flow is flrst attracted by the RFIM zero temperature fixed 
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point and approaches the standard first-order discontinuity fixed point for very large lengthscales only, see central 
and right panels of Fig 6. 

Close to Th the average value and the variance of the couplings start to grow as £^ . The field instead is characterized 
by a variance of the order £^ and an average (the configurational entropy) that is much smaller but grows faster, as 
(cif (T) — c)i^ (we use the standard notation for the RFIM exponents^). The point to set length is determined as 
discussed previously as the value of £ such that the typical renormalized couplings become of the order of the typical 
fields. This comparison leads to 



{ck{T) - c)2/(rf+^-2r,) 



where we have used the scaling relation x — 9 = ^((i + Ty — 277). By the same argument used in the previous section 
we also finds a diverging relaxation time-scale: 

pe 1 

t p c 1 

log r oc — ^ cx 



T (ckCT) - c)2e/('i+V2-2^) 

Since in d=3 the exponent 9 is very close to 1.5 and 77 ~ 27; (we indeed find by MK-RG x ~ 3), the relaxation time 
is expected to grow in this regime almost according to the so called Vogel-Fulchcr-Tamman law but with respect to 
the variable ck{T) — c instead of T — Tk, i.e logr oc [ck{T) — c)~^ . 

This critical behavior described above is valid until Ips becomes of the same order of ^ = {Th — T)^'^ {v is the 
RFIM length-critical exponent), i.e. if one does not approach ck close enough. Afterwards, for i!ps ^ ? we recover 
the behavior described in the previous section, which is due to the discontinuity fixed point. In order to understand 
the interplay between these two regimes we use scaling theory and assume that the averages of the coupling and of 
the field increase respectively as: 

This is corroborated by the numerical analysis of the RG flowSl. Since the scaling with £ of the average coupling and 
of the average field is and £''■ for the discontinuity fixed point, the behavior of the scaling functions has to be 
fi{y) oc yd--'^-^ and /2(j/) oc y''-^^ for large values of the argument y. This implies that eventually, very close to ck{T) 
for any temperature T < T^, one recovers the scaling discussed in the previous section but with singular pre- factors. 
In particular, the average coupling and the average field increase as ^s-d+i£d-i ^^^^ ^x-d^d^ leading to divergences of 
time and length-scales such as: 

£pS (X —- r , lOgT OC 



{ck{T)-c) ' ^ {ck{T) ^ c)(d-i) 

where v = v{d — 9 — \) + [d — l)(a; —1 — 6') ~ 2.125 in three dimensions (to obtain this result we used v ~ 2.25 
obtained from previous Migdal Kadanof RG computations^!^, using instead the numerical resullj^ v = 1.37 we have 
V ~ 1.69). Note that close to Th the MCT relaxation time also diverges with a generalized Vogel-Fulcher law^, thus 
a quite intricated set of cross-overs takes places close to the endpoint of the RPGT line. 

The value of the discontinuous jump of the overlap at the transition scales as the magnetization jump for the RFIM, 
i.e. 

A9 oc T"'' = [Th - Tf 

We recall that because of the scaling relation valid for the RFIM j3/v = d — x^. Since in three dimension^ (3 ~ 0.017 
the transition is expected to still look first-order like even close to Th, except if one is able to reach temperatures 
extremely close to Th (and for very large system sizes). The cross-over between the two regimes takes place for 
£ps ~ ^, which corresponds to concentrations of pinned particles such that {ck{T) — c) oc {Th — TY'^^~^^ (where in 
three dimensions v{x — 0) ~ 3.375 using u ~ 2.25, or ^{x — (?) ~ 2 using 1/ ~ 1.37). 



V. MAPPING TO THE RFIM AND SIMULATIONS OF THE THREE DIMENSIONAL ISING MODEL 

WITH RANDOM PINNING 

We have shown in the previous section that the MK-RG equations for the replica field theory map. after some 
approximations, to the ones valid at zero temperature for a peculiar random field Ising model characterized by a field 
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FIG. 8: Overlap with the reference configuration as a function of time for the 3D Ising model at T = 4, /i = 0.12 on 
a cubic lattice with iV = 8 * 10'^ spins (J is set equal to one). Increasing the fraction of frozen spins a much larger 
time is needed to escape from the high overlap phase. Inset: overlap distribution slightly above and below c* (T = 4, 

h = 0.09, N = lO^). 



that is negative and uniform but in randomly distributed sites, where is very large and polarizes the spins in the up 
state. The RPGT discussed previously corresponds in this setting to the transition from the negatively magnetized to 
the positively magnetized state, which becomes stable for a large enough fraction of spins pinned in the up position. 
When we established the mapping we neglected fluctuations around the saddle point and, as a consequence, we 
obtained the zero temperature MK-RG eqs for the RFIM. A heuristic, simplistic but suggestive way of taking into 
account the role of (at least some) fluctuations and to study the dynamical evolution for the pinned system^ is to 
consider the same spin model but at finite temperature. Just after pinning, and below ck{T), the liquid has a high 
overlap with the initial equilibrium configuration; it remains in this state for a long time (diverging for c — >■ ck{T)), 
but it eventually decorrelates from the initial condition and equilibrates toward the low overlap state. Similarly, 
consider an Ising model below and take an equilibrium configuration in presence of a positive magnetic field h] 
then freeze a fraction c of spins at random and revert the field to —h. Now the original configuration is metastable 
for c smaller than a critical value c*, and stable for larger ones. Despite some important differences, the situation is 
very reminiscent of the one we discussed before for supercooled liquids with h playing the role of the configurational 
entropy, the positive and negative magnetized states playing the role of high and low overlap ones. The life-time of 
the positively magnetized metastable state, which diverges for c* — c — > 0, is analogous, mutatis mutandis, to the 
relaxation time that takes to the liquid to decorrelate from the initial configuration and that diverges for ck — c ^ 0. 
As benchmark for the super-cooled liquid case, we have therefore studied the behavior of the overlap with the initial 
positively magnetized configuration as a function of time in the corresponding Ising model, see Fig. [H We indeed 
found that by increasing the value of c the system spends more and more time in the metastable high-overlap state 
before equilibrating at low overlap. The timescale for relaxation at low overlap diverges at a critical value c*{T,h) 
with the same law discussed before in the glass transition case: logr oc l/(c — c*)'^~^. Actually, the curves in Fig. [8] 
reproduce qualitatively well^ the ones obtained in super-cooled liquids by Berthier and Kob^S. This is an indication 
that the RFOT theory, and the related physical mechanism we have unveiled in and in this work, provide a 
viable explanation for the mysterious and generic dramatic slowing down of the dynamics, which has been observed 
in simulation when pinning particles of super-cooled liquid a^^'^^i^^'^^ . (The major important difference between this 
simple Ising model and real liquids is that the original reference configuration is not an equilibrium one. The main 
consequence is that the system is not automatically at equilibrium after having frozen out particles; thus, relaxation 
and equilibration time are not the same, contrary to the case of super-cooled liquids.) In the Ising model case we 
can also study the equilibrium of the system in proximity of the critical fraction of blocked spins c* . As it is usual 
for a first order transition, we find a bimodal distribution for the order parameter, which is the overlap between the 
configuration used to block the spins and any equilibrium configuration, see inset of Fig. [S] The transition point 
corresponds to perfect bi-modality, whereas instead close to the transition one of the two peaks eventually disappears 
for large enough system size. This, as we shall discuss in more detail in Sec. 7, suggests a way to study thoroughly 
RPGTs: exactly at ck{T) the average probability distribution function of the overlap with the reference configuration 
used to pin the particles is expected to show a double peak structure of the same kind, while a single peak structure 
with a peak at low (high) overlap should characterize the probability distribution function at c < ck(T) (c > ck{T)). 
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VI. COMPARISON WITH OTHER APPROACHES 

In this work and in Q we considered one more control parameter besides the temperature in order to favor the 
glass phase and in order to give rise to a line of phase transitions having the usual glass transition as one of the two 
endpoints. The very same idea is also at the root of two other precedent approaches: the ones developed by Franz 
and ParisiiS and by Chandler and Garrahan^. Despite their resulting phase diagrams look similar at first sight, there 
are actually very important differences, as we explain in the following. 

A. Space-time phase transitions and the s-ensemble approach 

Chandler, Garrahan and co-workers proposed to enlarge the number of control parameters for a super-cooled 
liquid by biasing the so-called activity of a given dynamical historyii*^i^. The idea is that the interesting dynamical 
features in proximity of the glass transition, like the local inhomogeneities in space and time called dynamic 
heterogeneities, show up in the space-time configurations as a coexistence phenomenon between regions showing high 
and very low local activity. A Lcgendrc transform of the action with respect to an external variable, s, coupled to 
the activity, is able to provide a quantitative analysis of this phenomenon and to actually reveal the occurrence of a 
first-order transition in the s,T plane. The resulting phase diagram is similar to ours in Fig. [1] besides the (relatively 
important) difference that the glass transition of ordinary unbiased super-cooled liquids takes place at zero, and not 
finite, temperature^. Even though the two approaches share this and other similarities, they differ in the detailed 
predictions and the physical content. 

The most important difference is that the transition on the line in the s,T plane is first-order (in space-time): the 
dynamical behavior changes discontinuously across the transition without showing any divergence of the relaxation 
time. Approaching the glass transition of unbiased super-cooled liquids by decreasing the temperature is therefore 
very different from approaching the transition line by increasing s. In our approach, instead, the RPGTs have 
the same critical properties for c t of the usual glass transition obtained by decreasing the temperature. In 
particular the relaxation time diverges approaching the RPGTs by increasing c. Other substantial differences can 
be found studying limit cases and dimensionality dependence. For instance, it has been shown that the space-time 
phase transition is present also in one dimensional models^, whereas as we pointed out previously the lower critical 
dimension for the RPGTs is two. 



B. The e-coupling approach: a detailed comparison 

An approach to the glass transition more similar to the procedure of pinning particle is the one developed by Franz 
and Parisii^. It consists in introducing a bias, e, in the thermodynamics of a glassy system to favor configurations 
correlated with a reference one, chosen from the equilibrium measure at temperature T. Both the e-coupling and the 
pinning procedure are equally aimed at reveal the complex multi-state structure in the phase space of glassy systems. 
In the e-coupling case, as for RPGTs, for every temperature in a finite range Tr- < T < T* it exists is a critical value 
es(T) for which a thermodynamic transition occurs in the biased system. The resulting phase diagram in the e,T 
plane is very much reminiscent of the c, T one (compare the phase diagram in (lo| with Fig.[T]). This notwithstanding, 
the differences between the two approaches are several. The most important one is that the transition at (T) is not 
a glass transition but is first-order with non-zero latent-hea1>i^ and, hence, is not accompanied by the divergence of 
the equilibrium correlation time. 

The physical reason behind these differences is that in the pinning case the metastable state correlated with the 
reference configuration is just one out of the very many in which the system can freeze. It does not differ from them 
neither for its internal energy, entropy nor for other properties. Instead, the e-coupling, introduced by an external 
term in the Hamiltonian, makes the state correlated with the preferred configuration very special and different from 
all the other ones. For example, the overlap among configurations belonging to the favored minimum is different from 
the overlap among configurations in any of the other available minima. The internal entropy and even the energy 
of the favored minimum become extensively different from the ones of other metastable states. Thus, whereas by 
increasing c one reduces the number of possible metastable states until the RPGT where the configurational entropy 
vanishes, this phenomenon is pre-empted in the e-coupling case by a first-order transition between the preferred state 
and all the other possible ones. A recent numerical study has shown that by constraining the overlap to a value 
intermediate between qea and zero, the system shows spinodal decomposition, as expected for a first order phase 
transition^. This is an indication that the phase transition predicted with the e-coupling indeed takes place (at least 
for a finite value of e) . 
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FIG. 9: The Franz and Parisi phase diagramiS (e^, esp,^ and es_p,^ hnes) for the p — 3 spin model has been 
completed by the introduction of the line of glass transitions for the metastable liquid phase (ck) and the Mode 
Coupling transition line (ed)- The inset shows a zoom of the low temperature region of the phase diagram where 
four ideal experimental paths {A, B, C and D) are superimposed to show the possible variety of transitions that it is 

possible to encounter using the e-coupling set-up. 



In order to analyze the difference between the two approaches we computed the coupling ek (T) that would be able 
to induce an entropy crisis of the liquid phase. It would represent, in the e-coupling framework, the analogous of 
the glass transition line ck{T). Starting from the results of [loj for the spherical p{= 3)-spin model and imposing 
that the low overlap phase must have zero configurational entropy, we obtained the full black line in Fig. |9] (See 
Appendix C for the computation). This glass transition line entirely lies in the region where the liquid phase is 
metastable. In consequence, it is irrelevant as far as equilibrium dynamics (or thermodynamics) is concerned, except 
for the initial point {exiTK) = 0,Tx) where it coincides with the glass transition in temperature. In the phase 
diagram of 1 1011 there are also two more lines very much reminiscent of the dynamic transition lines Cd{T) and c^^ (T) 
obtained in |41| and in this paper, and reported in Fig.[TJ However, also in this case appearances are deceiving. In the 
e-coupling phase diagram, only the line esp,^ reported in Fig.[9l representing the spinodal points of the liquid phase, 
is probably the analogous of Cdf below the e^ line^S. At variance, the line es p^^ is the line of spinodal points of the 
reference amorphous metastable state. These points are in general different from the points where generic amorphous 
metastable states appear in the liquid phase, which correspond to the usual dynamical transition a la Mode Coupling, 
i.e. the Cd{T) transition line represented in Fig. [TJ We have explicitly computed the MCT line, ed{T), and plotted 
it in Fig. [5] (It was obtained by finding the smallest e such that equilibrium of the liquid phase is characterized by 
the presence of a non-trivial solution with qi ^ qq and m = 1, in case of a 1-RSB ansatz of the replica matrix.) The 
distinction between the two approaches is now evident both in the statics (the transition is first-order and not RFOT 
like) and in the dynamics (the region where aging is expected, i.e. the region below ed{T), is much smaller in the 
e-coupling case). 

The phase diagram for the e-coupling reveals a quite intricate spectrum of transitions when e increases (see inset of 
Fig. [5]). When e rises at high temperature (line A), only a first order phase transition occurs at es{T) between a 
simple liquid phase and the amorphous state chosen by the reference configuration. This transition is preceded and 
followed by the spinodal points of the two phases. For temperatures larger but almost equal to the Mode Coupling 
temperature of the free system, two different cases are possible: at higher temperature (line B), increasing e, a well 
defined amorphous reference metastable state forms before the occurrence of a first order transition between the liquid 
phase and the amorphous state. However, in this case, a further increasing the coupling between replicas produces 
in the liquid phase (now metastable) first the formation of a large number of metastable states (the ergodicity of 
the liquid phase breaks) and finally the disappearance of these new states and of the liquid phase itself. At lower 
temperatures still larger than Td (line C), the first order transition starts to occur after the ergodicity breaking 
transition of the liquid phase. This implies that eventually a glassy behavior in the equilibrium dynamics becomes 
visible. At even lower temperature (line D), finally below Td, the liquid phase is non-ergodic for the free system, 
yet the thermodynamic transition occurring when e increases continues to be a first-order transition. Beyond this 
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transition, following the metastablc liquid branch, it is now possible to approach at last the true thermodynamic glass 
transition associated to the vanishing of the configurational entropy. 



VII. ON THE DIFFERENCE BETWEEN RANDOM PINNING PROCEDURES AND ON POSSIBLE 

NUMERICAL AND EXPERIMENTAL STUDIES 

The procedure of pinning particles has emerged as a way to probe the complex structure of the phase space of 
supercooled liquids and to study the occurrence in real systems of ideal glass transitions at ck{T). It can be used to 
compare and contrast the existing dynamical or thermodynamic theories of glass transition. In this discussion section 
we address several general issues related to it. 

A. Pinning from an equilibrium configuration versus pinning from a high temperature (random) 

configuration 

A natural question is how crucial is that the configuration used to pin particles is an equilibrium one. Actually, 
pinning from a high temperature (random) configuration is a procedure that has been already studied in the past 
in connection with the physics of glass-forming liquids in porous materials, both numerically in simulations^ and 
analytically within MCT^. 

Our study, based on mean-field glassy systems, show that the two procedures lead to a very different physical behavior. 
In the case of RPGTs the phase diagram displays a transition line with an endpoint at a high temperature Th- This 
implies that it is possible to follow a continuous path in the phase diagram connecting the liquid and glass phases 
without crossing any phase transition, see Fig. [TJ The main reason that makes this (at first sight) surprisingly fact 
possible is that the ideal glass phase is "simple": it just consists in the equilibrium configuration from which particles 
are pinned and vibrations around it. This clarifies why it is not paradoxical that liquid and ideal glass are connected 
by a continuous path and why the ideal glass phase can be sampled easily. In the case of pinning from a random 
configuration the ideal glass phase is instead "non-simple": there is no way to sample it easily and it is separated 
from the liquid phase by a transition line. This is expected on general grounds since this way of pinning particles 
(or spins) introduces quenched disorder that is not correlated to any particular state and looks completely random 
(besides possible short-range correlations). Thus, it does not favor any state; instead, it introduces more randomness. 
The phase diagram, see Fig. [2l contains a transition line that first, at low c, consists in RFOT-like transitions and 
that is actually similar to the corresponding one of RPGTs (except for the differences in the low-temperature phase). 
However, after the critical point, the transitions change nature and become continuous. This is interesting because 
it predicts that in this case, except at very low pinning, the physical behavior and the criticality are quite different 
from the ones displayed or expected for unconstrained super-cooled liquids. In particular there is a very large regime 
in which the transition is continuous. Instead of observing the characteristic two-steps behavior for the collective 
correlation functions one should instead find an overlap qi that emerges continuously from the long-time limit value 
go- Moreover, approaching the critical point the spin-glass susceptibility diverges and the relaxation time should 
increase in a power law way with the distance from the transition, as it happens for spin-glasses. These are quite 
remarkable and unusual features for super-cooled liquids. Actually, in a very recent work Karmakar and Parisi^ 
studied the behavior of liquids in presence of particles pinned from a high temperature (hence, essentially random) 
configuration and found evidences of a transition different from RFOT and possibly continuous. It is certainly worth 
investigating numerically the properties of this transition and determining its universality class. Along the same lines 
a detailed mean-field study of the general case in which particles are pinned from a configuration equilibrated at 
temperature T' different from the reference temperature T would be useful. 

Finally, we would like to mention that the difference we obtained between the physical behaviors induced by the 
two pinning procedures was in contrast with a prediction obtained by Krakoviack^ through an extention of MCT to 
pinned particles. This was very surprisingly since MCT is thought by many (and by us) to be an essential part of 
RFOT physics. Thus the predictions we found by mean-field, phenomenological and RG arguments were expected to 
be in agreement with a full fledged MCT computation. In a recent work Flenner and Szamei^ developed a different 
MCT theory for pinned particles in which static and dynamics are treated on equal footing. Their analysis is in 
complete agreement with ours, thus solving, at least partially, the difficulty mentioned before. 
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B. Looking for RPGTs in real systems: physical phenomena and finite size scaling 

All important open issue of the thermodynamic (RFOT) approach to the glass transition is to ascertain the existence 
of a true transition and determine its properties. This seemed an hopeless task: because of the dramatic growth of 
the relaxation time and the modest increase of the correlation length studying the critical regime is very difficult 
and, since the ideal glass phase cannot be reached at equilibrium, showing evidences of an ideal glass transition very 
challenging. The study of RPGTs provides a new and promising way to tackle this problem and circumvents the 
difficulties cited above. 

The main idea is that analyzing and proving the existence of the upper part of the RPGT line is doable (even 
if numerically challenging) by numerical simulations and, possibly, experiments with optical tweezers. From the 
knowledge gained in these studies and by extrapolation to lower temperatures one can obtain new information on the 
putative glass transition taking place at Tk for unconstrained super-cooled liquids. 

In the following wc discuss peculiar physical phenomena induced by RPGTs and propose ways to test the existence 
and the properties of these transitions. By approaching ck{T) the point-to-set length of the pinned system diverges 
as a power law, while the relaxation time is conjectured to experience a super-Arrhenius behavior. Thus, the first 
basic test is that by increasing c one has to find a really dramatic increase of the relaxation time, analogous to the one 
obtained decreasing the temperature for usual super-cooled liquids. This has been actually observed in very recent 
numerical simulation o^°i^^ . However, this is evidently not enough to assess the existence of a transition at any ck{T). 
The great advantage of RPGTs is that the glass phase is "simple" : it is characterized by the amorphous order imposed 
to the system by the configuration used to pin particles. Thanks to this property, probing the equilibrium regime 
beyond the transition is in principle an easy task since it is straightforward to start from an equilibrium configuration 
(the initial configuration used to pin particles). Moreover the equilibrium decorrclation time is not related to collective 
rearrangements and is given by a microscopic time-scale, which remains finite when ck{T) is approached from above. 
At variance, as recalled above, approaching the glass transition from below (using concentrations of pinned particles 
smaller than ck{T)) gives the usual super-Arrhenius divergence of relaxation time. Hence, a clear dynamical signature 
of the transition, and a way to pinpoint its location in the phase diagram is observing a non-monotonic behavior of 
the equilibrium decorrelation time when c increases and crosses ck (T) . 

Dynamical observations are not enough to certify the existence of RPGTs because, obviously, directly observing the 
occurrence of a divergence is impossible and the non-monotonic behavior of the relaxation time could be the indication 
of a simple crossover or an avoided transition. Moreover, in real particle systems, the microscopic time can be quite 
large because by pinning particles one freezes a lot of vibrations that are instrumental in providing fluctuating free 
volume that helps single particle motion. This effect can mask the presence of a non-monotonous behavior of the 
relaxation time (a possible solution to avoid this spurious and quite annoying byproduct is, instead of pinning, to 
restrain the motion of a fraction of particles by letting them only vibrate around fixed positions^. 
A crucial extra piece of information can be obtained by a systematic study of the thermodynamics and in particular 
of the overlap q between typical configurations extracted from the equilibrium measure. Indeed, the presence of a 
RPGT would imply a sudden change in the distribution of the overlap exactly at Ck{T) from a distribution with a 
low overlap peak, characteristic of the liquid phase where equilibrium configurations are typically far apart one from 
the others, to a distribution with a high overlap peak, characteristic of the trivial glass phase correlated with the 
initial equilibrium configuration. This sudden change is smeared out in systems that can be numerically studied in 
practice because of their finite size. This actually is not necessarily a limitation; a thorough study of finite size effects 
can provide very valuable information on the critical properties as it happens for standard phase transitions. 
We now present our predictions for finite size effect based on the RG study presented in this paper and in 
we expect three different regimes, which are controlled by the distance from critical point Cf^Th belonging to the 
universality class of the RFIM. Let consider the situation in which temperature is close to Th, the case relevant for 
numerical simulations. Both below and above Th, increasing the system size L for L < ^ (with ^ the RFIM correlation 
length), the scaling theory of the RFIM implies that the distribution of the overlap is characterized by two peaks 
whose distance scales as The bi-modality is due to disorder fluctuations and to the simultaneous presence 

of positively and negatively magnetized regions. Each single peak is broadened by thermal fluctuations on a scale 
l/l{d+n-2)/2 ^ i/Ld--^, Exactly on the transition line (below T,,) or on its continuation (above T^) both peaks 
shrink and increase with i in a similar way. Instead, if the system is beyond or before the transition line one, of 
the two peaks get depressed and eventually disappears for large is. Because of the value of the exponents in three 
dimensions^, d — .t ~ 0.011 and {d + rj — 2)/2 ~ 0.75, the positions of the peaks decreases very slowly with L. In 
consequence, except if one is able to reach very large system sizes, the variation of the peak positions with L should 
be almost invisible in practice. When L becomes of the order of ^ the behavior starts to change depending whether 
T > Th or T < Th. In the former case the distribution ceases to be bimodal, hence signaling the absence of the 
transition. In the latter, instead, the behavior crosses over to the one of long-range ordered systems, i.e. the distance 
between the two peaks becomes L-indcpendent, of the order of 1/^''"'^, and their width decreases as l/L'^/^. 
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The concomitance of the dynamic and thermodynamic features discussed above and detailed finite size scaling studies 
of the probability distribution functions of the overlap provide a new way to ascertain (or disprove) the existence of 
RPGTs for any Tk < T < Tji and determine its position ck{T) in the phase diagram. 

Another interesting and promising research direction that we do not cover here since it goes beyond the scope of this 
work is studying aging in presence of random pinning. We refer the interested reader to our paper— for a mean-field 
analysis. 

C. Dimension and protocol dependence 

Using the insights gained by our RG analysis one can obtain interesting predictions concerning the dimensionality 
and the protocol dependence. These were already discussed in Q. We report them here for completeness. 
Because of the mapping on the RFIM, our RG procedure predicts no transition, just a cross-over, in two dimensions. 
Thus, the lower critical dimension for RPGTs is two. Before embarking in numerical tests of this prediction it would 
be interesting to understand the sharpness of this cross-over; we leave this issue for future studies. 
Of the other pinning procedures introduced in the literature, two- walls, one- wall and the cavity ones, only the first can 
lead to a phase transition akin to the RPGT studied in this work. The reason in the cavity case is clear, the system 
has a finite size, for the one- wall case, see [59| . It is important to realize that in the two- wall case the free system 
has effectively one dimension less, e.g. in three dimensions the sandwich between the two walls is two-dimensional 
(the system is only infinite in the two dimensions parallel to the wall). Because of this, the lower critical dimension is 
shifted up of one unity and no transition is expected in three (and lower) dimensions. This geometry has been studied 
very recently and, indeed, the corresponding results are compatible with this predictior^. 

VIII. CONCLUSION 

In this work we have presented a detailed study of Random Pinning Glass Transitions that includes mean-field 
and RG analysis. We have obtained the critical properties and discussed possible numerical and experimental tests. 
Random pinning appears to provide a new angle to tackle the glass transition problem. We are optimistic that in the 
future, by progressively investigating the physical behavior and the finite size effects for randomly pinned fluids and 
studying different ways of pinning particles, a lot of substantial progress will be made. 
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Appendix A: Entropy vanishing transition for a pinned particle system 

One of the advantage of pinning particles from an equilibrium configuration in a supercooled liquid system is the 
possibility of inducing an entropy vanishing transition at temperature higher than Tk. The existence this kind of 
transition in real systems and its properties is matter of a longstanding and still open debate. In this appendix we 
present a short argument that shows in more details why the transition induced by pinning particles is necessarily 
associated to the vanishing of configurational entropy. 

The partition function for a liquid replicated m times and where a fraction c of particles are all constrained to be in 
the same replica can be schematically written as follows: 

Z'" = ^exp(-Af/3/,)exp(-7V/3(m - l)/.es(/a)) = exp(iV5(c, T)) (Al) 

a 

where the sum runs over the metastable states indexed by a. Only the first replica is allowed to freely explore these 
metastable states. The other to — 1 replicas are characterized by a reduced free-energy /res(/) and configurational 
entropy Sresif) due to the degrees of freedom of the remaining free particles. According to the usual procedure, we 
can rewrite the replicated partition function as an integral over the levels of free-energy density. Doing this we have 
to consider two cases, the first one is the case where the free degrees of freedom of the constrained replicas are still 
free to explore an exponentially large bunch of metastable states, in this case we have 

= / dfcM-N[3f + Ns,{f)) exp(-7V/3(TO - l)r"^(/) -I- N{m - l)sr'(/)) (A2) 
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which in the thermodynamic hmit, using to the saddle point method, becomes 

= exp(-7V/3r + Ns^iD) exp(-7V/3(m - l)r^^(r ) + iV(m - 1)^^ (/*)) (A3) 

where /* is such that 



df 



(A4) 

/* 



when m — > 1. In the second case, the constrained rephcas are stuck in the amorphous configuration given by the 
pinned particles, hence the partition function simply reads 

d/exp(-7V/3/ + se(/)) exp(-iV/3(m - l)r'{f)) (A5) 

and in the thermodynamics limit it becomes 

= exp(-iV/3/* + ,s,(/*)) exp(-iV/3(m - (A6) 



where (jA4p also holds. The transition between the liquid and the amorphous crystal happens when the total free- 
energy associated to these two partition functions are equal. This automatically gives at the transition the condition 
3'^"%/*) = 0. 



Appendix B: Inducing a glass transition in the e-coupling framework 

The e, T phase diagrant^^ obtained studying the thermodynamics of a spherical p-spin system in presence of an 
e-coupling between equilibrium configurations apparently shares some features with the c, T phase diagram obtained 
by pinning particles from an equilibrium configuration. In this appendix we re-propose the derivations of the two 
spinodal lines and the thermodynamic transition line of this phase diagram. These transitions are actually qualitatively 
different from the Mode Coupling transition line and the glass transition line obtained by pinning particles from an 
equilibrium configuration. Then, for an actual comparison with the c, T phase diagram obtained by pinning particles, 
we propose an original derivation of the analogous of the Mode Coupling transition line and of the glass transition 
line in the e, T phase diagram. The novel e, T phase diagram comes out to be quite different from the c, T one. A 
detailed comparison between the two final phase diagrams is in the discussion section of this paper. 
The e-coupling procedure consists in biasing the thermodynamics of a spin system by an external field pointing 
towards an equilibrium configuration C^e/ at temperature T. We restrict here the discussion to the case where the 
temperature of the reference configuration Cref is equal to the temperature of the coupled system. If aref are the 
spins in the positions given by the reference configuration, the partition function of the biased system is 

Z[e, T, aref] = cxp(-/3i?jM + /3et7a„/) • (Bl) 
Averaging over the possible aref and the quenched disorder J it is possible to obtain the free-energy of the system 



^t'' ^] = Z[T] ^ eM~PHj[aref]) log(Z[e, T, aref]) (B2) 

where 

Z[r] = ^exp(-/3i/j[a„/]) , 

and its Legendre transform 

V[q;T]=minT[e,T]+eq^ F[T] (B3) 

e 

where 

F[T] = -ilog(Z[T]) . 
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The last function V[q; T] represents the thermodynamic potential associated to every possible value of the similarity 
q between the reference equilibrium configuration and its coupled configuration. By definition we have imposed 
1^[0,T] = 0. When a finite e-coupling is imposed, the equilibrium of the system is given by the minimization of 
V^[(3';T] — eg, hence by the condition 

^ = (B4) 
dq 

The averages in (|B2[) are performed through the replica trick by the introduction of two series of n and m replicas of, 
respectively, the reference configuration and the coupled system. Suitable overlap matrices among replicas of the two 
series are consequently used within the following restrictions: a replica symmetric ansatz for the overlap matrix among 
replicas of the reference configuration, which means a restriction to the T > Tk region of the phase diagram, and a 
IRSB ansatz with parameter qi, go and m G [0, 1] for the overlap matrix among replicas of the coupled configurations. 
Finally the overlap between a specific reference configuration and every coupled configuration is settled to q, while 
it is assumed to be zero for all the other reference configurations. Under these assumptions the potential 1^[(7;T] 
corresponds to the following expression: 

V[q;T]= max ^-L U 13' f{q) ~ [{1 - m)f{q^) + mf{qo)] + "^^logil - q,) 

+ - log(f - (f - m)gi - mqo) + 7-^^ 1 (B5) 

m 1 — (1 — m)qi — mqo J 

where f{q) — ^q^, and the equations for the stationary points read as follows: 

dV 

7^=0: qo-q' = f'iqa)[l-il-7n)qi-mq„f (B6) 

oqo 



oqi 
am 



+ - q,) - ^/'(qo)[l - (1 - m)q, - mq^] = (B8) 



where f'{q) is the total derivative of f{q) with respect of q. 

The solution of this set of equations for every value of q gives the parameters qi , go and m to be inserted into 
(|B5[) for finally have the potential V[(7;T]. IniS it was found that this potential has a stable minimum in g = at 
high temperature and develops at temperature Td a second minimum at finite overlap qsA- Further decreasing the 
temperature, the value of the potential in the secondary minimum decreases and reaches zero for T = — Tk- This 
indicate the occurrence of a thermodynamic transition associated to the vanishing of the configurational entropy. 
The same features can observed for the system biased by the e-coupling term. In particular at fixed temperature 
increasing e, the new potential y[g;T] — eg has a single minimum at overlap gi(e,r) for < e < esp^^{T). At 
egp^^{T) a new minimum in correspondence of qni^jT) forms and remains metastable for esp^^{T) < e < es(r). 
At €s{T) a thermodynamic transition occurs, beyond which the equilibrium is given by the high overlap solution 
between the reference configurations and the configurations coupled to them. Finally in the large coupling regime, for 
t~esp {T) the minimum at low overlap q^ie^T) disappears. The €s(T) line is hence defined by (|B4[) which has two 
solutions: gL(e,T) and ^^/(e, T), and the condition {V[q-T]-eq)\^^ = {V[q]T] - tq\^. The esP,jT) and egp,^ (T) 
lines correspond to the e such that both (|B4p and its derivative (fiVlq^T]/ dq' = are satisfied and for each T this 
happens for two distinct values of e. 

In general the transition at e^ has no reasons to takes place where the configurational entropy vanishes for the coupled 
system. As we will see in the following, the single point where the two transitions coincide is e = 0,T = Tk, 
anywhere else the transition at Cg is a standard first order transition. To understand this issue, we have to look at 
the solution gi,go and m of the set of equations (|B6IB7IB8|) for the stationary point of T] in correspondence 
of the stable gi,(es,r). These parameters give information on the equilibrium characteristics of the ensemble of 
configurations coupled by e to the reference configuration in such a way that the overlap between them and the 
reference configuration is q. In general, for every chosen g, the parameters gi,go and m change. For very small g, 
the first two equations of this set have not a solution with gi ^ go for m = 1, the system coupled to the reference 
configuration is a simple paramagnet before the Mode Coupling temperature. When g increases and reaches q = q'^, 
the non-trivial solution gi ^ go for m = 1 appears ; this corresponds at the Mode Coupling transition. Beyond 
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this point the phase space of the coupled system is broken up into mctastable states with self-overlap qi and mutual 
overlap go- Finally, the point where the configurational entropy of the coupled system is zero corresponds to q = q* , 
the point where the solution of the full set of equations gives m = 1. Beyond this point the full set of equations for 
qi, qo and m gives the result qi ^ qo and to < 1, meaning that the IRSB phase has been reached. Hence the point q* 
is the extremal point of the IRSB region. 

For example in the 3-spin case, in correspondence of qL{^s{T),T), the system coupled to the reference configuration 
is very far from the limit q* of the IRSB region and, usually, even from the crgodicity breaking point q'^ of the 
Mode Coupling transition. For this reason the transition at typically represents a first order transition between 
a trivial paramagnetic phase to the amorphous state pointed by the reference configuration. To actually find the 
dynamic and the glass transition line id{T) and exiT) we need to find the e for which respectively the q'^ and the q* 
becomes the stationary point qL(e,T) of the potential F[g;T] — eq. In particular we found a line of eif(T) located 
beyond the line es{T) where qL{e,T) becomes metastable, meaning that the coupled paramagnetic phase is never 
stable when the entropy vanishing transition occurs. On the contrary, the dynamic transition €d{T) is found where 
the paramagnetic phase is stable only at the lowest temperatures and metastable at higher temperatures. In summary, 
in the thermodynamic stable phase of the spherical p-spin, the e-coupling procedure is never able to induce a glass 
transition associated to the vanishing of configurational entropy and only at very low temperatures it is able to induce 
the dynamic transition due to the crgodicity breaking of phase space of the coupled configurations. 



Appendix C: Entropy vanishing transition induced by pinned particles in the e-coupling apporach 

The overlap potential V[q] T] obtained in the e-coupling framework provides a valid approach to the glass transition 
problem. We apply here this description to the case of a p-spin model with pinned spins from an equilibrium configu- 
ration. We will pay particular attention to the amount of the configurational entropy exactly at the thermodynamic 
transition ck(T) to show that this transition is always an entropy vanishing transition. 

When spins are pinned from an equilibrium configuration the potential between replicas coupled to have overlap q has 
a structure similar to the one presented in (jB5[) with f(a) — i(c-f (1 — c)^)^: V[q;c,T] = max^^ ,jq.,„ V[gi, ™; 9; c, T] 
where 

V[qi,qo,m;q;c,T] = 

^ ^2/32/(g) - - + mfiqo)] + (f - c)^^ log(l - q^) 



2/3 

+ (1 - c)l log(l - (1 - m)gi - mqo) + (1 ~ c) ' f 1 . (CI) 

TO 1 — (1 — mjqi — mqo J 

The points ck{T) where the static transition occurs satisfy the condition V[qL] c, T] — V[qH] c, T], where (7l(cj T) and 
qH{c,T) are such that dV/dq\q^ ~ and dV/dq\qjj = 0. Note that, first, T^[(7;c, T] is the result of a maximization of 
the action with respect to qo, qi and to. The equations for finding the stationary point in gi, go and to have to be 
satisfied by the to = 1 condition when the glass transition occurs. Second, this condition is equivalent to impose that 
the configurational entropy Sc = dV/dm\m=i,q=qL of the coupled configurations vanishes at the transition. We now 
show that the point where Sc vanishes coincide with the point V[qL; c,T] = V[qH', c,T] where the static transition 
occurs. 

The stationary point conditions over qi , go and m with the condition m = 1 determine the glass transition point and 
read as follows: 

~ P'l [[c+{l - c)q^r c)qon (C2) 



(l-gi)(l-go) " 2 



A+{l-c)q,Y (C3) 



(l-9o)2 2 



^2 



^ [[c + (1 - c)gi]3 - [c + (I - c)go]^] + (1 - c) lo; 



-f(l-c)i^l^M(i^^|^l±^=0. (C4) 

(1 - (igY 

The solution gf go to this set of equations has to be found in correspondence of the stationary points gL,i/(c, T) 
that satisfy the following condition 



P^hc+[l-c)qL,H?^-^^ (C5) 
^ 1 - 9o ' 
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obtained by imposing dV /dq] 



1L,H 



0. 



We must show that the glass transition point corresponds to the point V[qi ^q^ ~ l;gL;c, T] — V[qi ^ q^ ,m = 
1; qu', c,T] where the potential shows a thermodynamic transition with a discontinuous jump in the overlap q from q^ 
to qn- 

Two tentative solutions for the set of equations (|C2p . ()C3|) and (|C5|) are qu = q^ = Qi and qL ~ q^ qi ■ Indeed, 
when 9o = 'Jl.Hj the conditions (jCSp and (jCSP correspond to each other. Moreover the relation (|C2[) is identically 
satisfied for the solution q = qn while it remains as a constraint to = qi{qL, c, T) for the solution q = q^: 



9f 



r, jxr, ^ =P\ [[c + (1 - ]^ - [c + (1 - c)g,.]2] . 



(C6) 



Finally, the condition for the static transition V[(j'f , , m = 1; q^; c, T] = V[gf^ , q/f jUi = 1; qn', c,T] is 



1 

2^ 



^[c + (1 - c)g^]3 + (1 - c) log(l - qn) + (1 - c)^^ 



1 

"2^ 



^[c + (1 - c)gL]3 + (1 - c) log(l - qL) + (1 - c)^^ 



(C7) 



On the other hand, for the low solution, q ^ qL, the configurational entropy reads 

^2 /1_.,L 



[[c + (1 - c)q^r - [c + (1 - c)^^]^] + (1 - c) log (^-7^) + 



(C8) 

Comparing the configurational entropy expression with (|C7p . it is evident that we would have Sc = at the static 
transition, if ~ qn- Indeed, this is exactly what happens: coming back to (|C6|) and using (jCSp . we obtain that 
qi obeys to the same equation (|C5|) as qn- Hence we showed that the static transition is always accompanied by the 
vanishing of the configurational entropy. 

This particular symmetry in the stationary conditions for the potential V[qi, qo, to; q; c, T] is obviously also present in 
the free p-spin model that can be obtained from the analyzed case putting c = 0. This is at the origin of the usual glass 
transition at Tr-. At variance, this symmetry gets immediately broken when an external e coupling between replicas 
is added to the Hamiltonian. In this case, the previous set of values for the overlap qn = q§ = q^ = qi ^ qL = qt 
does not constitute a solution for the problem. Instead slightly different values of qs characterize qn, qo and the 
high overlaps qi in correspondence of qL and qn- Moreover an explicit term —eq makes the condition V[qf ,(7q,to ~ 
1; qL', c, T] = V[q^ , q^ , m = 1; qn', c, T] unmistakably different from (|C4[) where e does not appear. This indicates that 
the static transition induced by the introduction of any finite coupling between replicas is different both from the 
glass transition at Tn and the glass transition induced at higher values of the temperature by pinning particles. 
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Actually, this assumption is not necessary to explain the physical behavior of super-cooled liquids; the configurational 
entropy could eventually not vanish, even though all the experimentally accessible physical behavior may be explained as if 
it were vanishing at a finite temperature Tk- 

Actually in this case there is another phase entering in the game; it occupies the right part of the phase diagram, the 
blue region in the bottom panel. We call it frozen liquid. Like the amorphous crystal phase, it is constituted by a stable 
minimum in the free-energy landscape and can be obtained from the liquid phase without crossing any thermodynamic 
singularity. On the other hand, it is separated from the glass phase by a continuous transition. Hence metastability and 
nucleation phenomena does not show up in this case contrary to the other one. Finally it is not characterized by any 
particular amorphous order but it is by selected by a completely random configuration, i. e. a liquid configuration at infinite 
temperature frozen by the pinning particle procedure. 

More precisely it is the microscopic length scale corresponding to the first peak in the radial distribution function. 
We only fix the constraint for the sites in the center of the hierarchical lattice, corresponding to n = 1. The reason is that 
the other sites of the lattice, entering in the RG equations for n > 1, actually correspond to renormalized regions and not 
microscopic ones. 

Note that with a slight abuse of notation we use the same symbol c we used before for the fraction of frozen particles. 
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The effect of this contribution it to block a fraction c of spins in the up state, corresponding to the high overlap position. 
In a finite dimensional lattice one can take for the magnetic field any positive value larger than J times the connectivity of 
the lattice. In the case of MK-RG any field larger than 2J is enough. 

A complete check of this scaling hypothesis is quite difficult numerically since it requires to follow the flow on extremely large 
length-scales and very close to the endpoint. All evidences we have from the numerics corroborate the scaling assumption. 
Besides the trivial difference that asymptotic low overlap value is negative. 

Very recently, it was shown that is also possible that the lower endpoint of the transition line takes place for a value of s 
larger than zero^. 

This is reasonable if we think that the point where Cdf occurs is the spinodal of amorphous metastable states^ and in 
particular of the threshold metastable states, i.e. the most numerous metastable states that are placed at relatively high 
free-energy. In practice, for c < c^^ the non-equilibrium dynamics remains trapped in these threshold states, while for c > Cdf 
equilibrium in the reference metastable state can be recovered in a finite time. Below the Mode Coupling transition line, e^, 
that we will discuss better in the following, the liquid phase is essentially constituted by an assembly of these states, hence 
spinodal of the liquid phase coincide with the spinodal of the metastable states. Actually this conclusion is straightforward 
if metastable states disappear all in once. Otherwise we could imagine that the spinodal of metastable states at different 
level of free-energy progressively disappears form high to low free-energy. Even in this case equilibrium will be recovered 
only when the whole set of metastable states have disappeared hence Cd^ conceptually coincides with eg p^^ . However the 
definition of the line Cdf in strictly refers to the disappearance of threshold states. It would be interesting to perform a 
more detailed analysis to see whether the Cd^ line moves when we study non-equilibrium dynamics from initial equilibrium 
configurations at several intermediate temperature T' > T and not only in the case T' = oo^. 



